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Executive  Summary 


This  documentation  is  the  final  technical  report  for  work  performed  under  AFOSR  grant 
F49620-96- 1-0445  .  The  title  of  this  project  was  "Atomistic  calculations  of  defects  in  YAG"  and 
the  Principal  Investigator  was  Ravi  Pandey,  Physics  Department,  Michigan  Technological 
University,  Houghton,  MI  49931.  This  was  12-months  research  program  to  use  atomistic 
modeling  techniques  to  study  defect  energetics  in  YAG.  During  the  grant's  active  period,  the 
following  report  is  submitted  for  publication  : 

(i)  Atomistic  Modeling  of  Defects  in  YAG, 

M.  Kukla  and  R.  Pandey 

submitted  for  publication,  J.  Am.  Ceram.  Soc.  1998. 


Native  point  defects  in  Y3AI5O12  (YAG)  were  studied  in  the  framework  of  the  pair- 
potential  approximation  coupled  with  the  shell  model  description  of  the  lattice  ions.  For  the 
perfect  lattice,  a  new  set  of  potential  parameters  were  obtained  which  reproduce  structure,  elastic 
and  dielectric  constants  of  YAG  very  well.  The  calculated  formation  energies  for  native  point 
defects  predict  that  the  antisite  disorder  is  preferred  over  the  Frenkel  and  Schottky-like  disorder 
in  YAG.  The  calculated  value  of  the  distortion  caused  by  the  antisite  Yai  are  in  excellent 
agreement  with  the  EXAFS  measurements.  In  non-stoichiometric  YAG,  the  calculated  reaction 
energies  indicate  that  excess  Y2O3  or  AI2O3  is  likely  to  be  accommodated  by  the  formation  of 
antisites  rather  than  vacancies  in  the  lattice. 


Personnel  Supported : 


V 


Faculty  :  Ravi  Pandey 

Graduate  Research  Student :  Maija  Kukla 

Visiting  Research  Fellow  :  Manuel  Redo 

Interactions/Transitions  : 

1.  Visit  to  Keran's  group  at  Wright  Laboratory,  Dayton,  July  1997. 


New  Discoveries,  inventions,  or  patent  disclosures. 


None. 


Honors/A  wards 


None. 


Contents  : 


Executive  Summary 

Atomistic  Calculations  of  Defects  in  YAG . 3 


I.  INTRODUCTION 


Yttrium  aluminum  garnet  (YAG)  is  an  important  material  whose  tech¬ 
nological  applications  range  from  lasers  to  propulsion  systems.  For  example, 
the  Nd3+  doped  YAG  is  a  well-known  laser  material  whereas  the  YAG:Al203 
composite  is  found  to  be  an  ideal  material  for  high-temperature  structural 
ceramic  applications.1  Over  the  years,  most  of  the  experimental  efforts  have 
been  focused  on  the  understanding  the  optical2-6  and  magnetic7-10  proper¬ 
ties  of  dopant  ions  in  the  YAG  lattice.  However,  relatively  less  attention  has 
been  paid  to  investigate  its  properties  with  an  aim  to  exploit  its  potential  in 
high-temperature  applications.11-13 

Diffusion  and  defect  chemistry  of  pure  and  doped  YAG  have  been  a  sub¬ 
ject  of  a  few  experimental  and  theoretical  studies.3,14-18  An  analysis  of  elec¬ 
troconductivity  measurements  and  diffusion  coefficients  on  a  series  of  the 
crystals  with  the  garnet  structure  has  established  that  YAG  is  an  ionic  con¬ 
ductor.  However,  it  has  been  suggested  that  the  conductivity  exhibits  mixed 
ionic-electronic  character.14  for  temperatures  greater  than  800° C.  A  simi¬ 
lar  conclusion  was  also  achieved  from  the  study  of  transport  properties  and 
defect  formation  of  Ca-  and  Mg-doped  YAG.16 

A  theoretical  study  based  on  the  shell  model  has  only  considered  the  pres¬ 
ence  of  vacancies  and  interstitials  in  the  lattice  predicting  the  dominance  of 
vacancies  over  interstitials.17  However,  spectroscopic  study  of  stoichiome¬ 
try  deviation  in  YAG  has  shown  a  possibility  of  cation  antisite  substitution 
in  the  garnet  crystal  lattice.19  This  has  been  confirmed  very  recently  by 
the  EXAFS  measurements20  in  Y2O3  rich  YAG  identifying  the  local  order 
around  the  yttrium  antisite  atoms.  Furthermore,  the  atomistic  study  of  the 
defects  in  Yttrium  Iron  Garnet  (YIG)  by  Donnerberg  and  Catlow21  has  also 
concluded  that  antisites,  not  vacancies,  will  dominate  the  intrinsic  disorder 
in  this  material. 

Due  to  a  close  similarity  between  YAG  and  YIG  materials,  it  is  therefore 
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expected  that  antisites  are  likely  to  play  an  important  role  in  YAG  and  should 
be  considered  in  any  experimental  or  theoretical  study  of  its  defect  chemistry. 

In  this  paper,  we  seek  to  perform  such  a  task  reporting  the  results  of  a 
theoretical  study  of  native  defects  in  YAG.  Our  approach  is  based  on  the 
pair  potential  and  shell  model  description  of  the  ionic  interactions  in  the 
crystalline  lattice  for  which  we  will  determine  a  new  set  of  potential  param¬ 
eters.  This  set  will  then  be  used  to  study  intrinsic  disorder  and  reactions 
describing  deviations  from  stoichiometry  in  YAG. 


II.  METHOD 


Partitioning  of  a  crystalline  lattice  into  two  regions  is  the  basis  of  the 
cluster  model  which  is  applied  very  often  in  calculations  of  defect  properties 
in  ionic  materials.  The  main  idea  is  to  separate  a  cluster  (referred  to  as 
region  1)  from  the  rest  of  the  lattice.  It  is  also  assumed  that  the  perturbation 
caused  by  the  defect  in  the  crystal  is  weak  beyond  the  cluster.  The  cluster  is 
generally  consisted  of  a  limited  region  containing  defect  and  several  spheres 
of  its  nearest  surroundings.  A  number  of  ions  in  the  cluster  varies  from  a 
problem  to  a  problem  and  depends  on  the  type  of  defects  in  the  system.  In 
calculations  the  cluster  is  treated  explicitly  within  a  selected  method  whereas 
the  reminder  of  the  lattice  (represented  by  the  region  2)  is  treated  in  an 
appropriate  way.  Thus,  the  crystal  lattice  energy  can  be  written  in  the 
form:22 

E  =  E\{x)  +  E2(x,  <f)  +  Ez(£),  (1) 

where  E  is  the  total  energy  of  the  crystal,  Ei(x)  is  the  cluster  energy,  E3(<f) 
is  the  energy  of  the  outer  region,  E2(x,  £)  is  the  energy  of  the  interaction 
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between  regions  1  and  2  depending  from  both  ion  coordinates  of  the  cluster  x 
and  the  rest  of  the  crystal  £.  The  exact  expression  for  the  crystal  energy  can 
be  written  as  a  series  of  interaction  terms  containing  n-body  components. 
The  interionic  potential  model  used  in  the  present  work  is  well  described 
elsewhere,  that  is  why  we  give  here  only  main  ideas  and  do  not  stop  on  the 
details.  The  method  assumes  that  the  total  crystal  energy  can  be  given  by  a 
sum  of  all  pairwise  interaction^  in  the  lattice.  In  other  words,  all  terms  but 
two-body  contribution  are  neglected  in  the  series.  Then  one  can  write 

E  =  ££(>, }),  (2) 

i>j 

where 

T  Qi  ’  Qj/rij-  (3) 

E(i,j )  is  the  energy  of  the  interaction  between  ions  i  and  j,  with  charges  qi 
and  qj  located  at  the  distance  rtJ .  The  4>ij{rij)  is  some  suitable  pair  potential 
describing  short-range  interaction  between  ions.  The  latter  is  a  result  of  atom 
repulsion  due  to  their  electron  shell  overlapping  and  dispersion  attraction. 
For  ionic  materials  this  is  generally  chosen  to  take  the  form  of  a  Buckingham 
potential:23 

<f>ij(rij)  =  Aij  exp (-rij/pij)  -  Cij/rfj  (4) 

Aij ,  pij  and  C,j  are  empirical  parameters,  that  are  usually  fitted  to  crystalline 
properties. 

Displacements  of  ions  in  inner  region  1  are  obtained  explicitly  by  mini¬ 
mization  of  the  total  energy.  For  outer  region  2,  displacements  of  ions  are 
determined  by  the  electric  field  due  to  the  effective  charge  of  the  defect  in 
the  framework  of  Mott-Littleton  approximation.24  For  cubic  materials  such 
as  YAG,  the  ionic  displacements  are  isotropic  and  are  given  by:25,26 

4»  (,-r  •  G-'  ■  q) 
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in  which  6x  is  a  vector  of  displacements  in  x,  y  or  z  directions,  Dx  is  the 
corresponding  electric  displacement,  G  is  the  second  derivative  of  the  non- 
Coulombic  interaction  matrix,  and  e  is  the  trace  of  the  dielectric  constant 
tensor. 

The  interionic  potential  model  is  significantly  improved  if  one  takes  into 
account  the  electronic  polarization  introduced  by  the  charged  defects  in  the 
lattice.  This  can  be  done  by  the  shell  model  of  Dick  and  Overhauser27  de¬ 
scribing  valence  electrons  as  a  massless  charged  shell,  that  is  coupled  with 
core  of  the  ion  by  isotropic  elastic  forces.  All  formulae  are  still  valid  if  one 
rewrites  them  substituting  shell  and  core  charges  and  coordinates  instead  of 
rigid-ion  ones.  Both  shell  and  core  charges  and  coupling  constants  of  a  shell 
and  a  core  of  the  given  ion  are  the  parameters  of  the  theory  and  are  generally 
empirically  derived. 


III.  RESULTS  AND  DISCUSSION 


Experimental  studies  of  garnets  started  long  ago  since  the  garnet  struc¬ 
ture  was  originally  solved  by  Menzer. 28,29  However,  there  are  only  a  few 
theoretical  investigations  focussing  on  defect  chemistry  in  garnets. 17,30,21,31 
The  garnet  structure  belongs  to  the  space  group  Ia3d  (0^°).  The  cations 
are  all  in  special  positions  labeled  as  a,c,d  with  no  positional  degrees  of 
freedom,  while  the  oxygen  atoms  are  placed  in  the  general  positions  96(A). 
Yttrium  occupies  dodecahedral  24(c)  positions  whereas  there  are  two  differ¬ 
ent  sites  for  aluminum  ions,  namely  octahedral  16(a)  and  tetrahedral  24(d) 
in  the  lattice  (see,  for  example,32).  This  leads  to  the  formula  of  Y3Al2Al30i2 
which  is  commonly  written  as  Y3Al50i2.  The  unit  cell  is  large  and  consists 
of  8  formula  units  (160  atoms).  For  calculations,  the  oxygen  positional  pa- 
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rameters  in  the  lattice  were  taken  from  neutron  diffraction  measurements 
by  Prince33  and  are  given  by  x=-0.029,  y=0.053  and  z=0.151.  The  X-ray 
diffraction  study  by  Euler  and  Bruce34  reports  the  cation-oxygen  distances 
in  YAG  to  be  1.94  A  and  1.76  A  for  Al3+(a)-02-  and  Al3+(d)-02-,  respec¬ 
tively.  Following  Geller,35  the  most  accurate  value  of  the  lattice  constant  for 
the  stoichiometric  Y3AI5O12  is  12.000i0.002  A.  Several  other  experimental 
studies  find  the  lattice  constant  to  be  slightly  higher  which  may  be  due  to 
the  presence  of  excess  yttrium  in  the  lattice.35,36 


A.  Interatomic  Potentials 


For  interatomic  potential  model,  formal  ionic  charges  are  assigned  to  the 
host-lattice  ions,  thereby  considering  them  as  Y3+,  Al3+  and  02-  in  YAG. 
The  short-range  interaction  terms  considered  in  the  model  describe  cation- 
oxygen  and  oxygen-oxygen  interactions.  The  short-range  interactions  be¬ 
tween  cations  are  ignored  as  they  generally  become  very  small  for  large  sepa¬ 
rations  which  are  2.99  A  for  d(Y-Al(d));  3.67  A  for  d(Y-Y)  and  d(Al(d)-Al(d)) 
and  5.19  A  for  Al(a)-Al(a)  in  YAG. 

Note  that  the  parameters  for  the  02--02--interaction  are  assumed  to 
be  transferable  among  a  series  of  oxides  and  are  taken  from  the  work  of 
Catlow  et  a/.22,37,38  They  were  obtained  from  Hartree-Fock  calculations  of 
the  interaction  of  two  negative  oxygen  ions.  Since  the  accurate  representation 
of  ionic  polarization  of  a  crystalline  lattice  is  known  to  be  very  important 
for  defect  modeling,  we  have  considered  all  of  the  constituent  ions  of  YAG 
(i.e.  Y3+,  Al3+  and  02'  ')  as  polarizable  ions  in  the  lattice.  Model  parameters 
are  then  fitted  to  experimentally  known  YAG  crystal  properties  such  as  the 
structure,  elastic  and  dielectric  constants  and  are  given  in  Table  1.  (For 
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details  of  the  fitting  procedure,  we  refer  to  the  reference39).  As  it  can  be  seen 
from  Table  2,  the  model  parameters  reproduce  the  crystalline  properties  of 
YAG  very  well  thus  providing  us  a  sound  basis  for  extending  the  model  to 
defect  calculations. 

Recently,  Bush  et  al.40  have  obtained  el  consistent  set  of  pair  potentials 
empirically  fitted  specifically  to  the  experimentally  measured  lattice  proper¬ 
ties  of  crystal  oxides  including  YAG.  This  set  is  in  addition  to  the  potential 
set  used  earlier  by  Schuh  et  al.17  which  was  fitted  to  AI2O3  and  Y2O3  struc¬ 
ture  properties.22,41  Although  these  two  sets  of  interatomic  potentials  yield 
good  agreement  with  the  experimental  values  of  the  lattice  structure  and 
elastic  constants  of  YAG,  they  do  not  yield  the  low-  and  high-frequency  di¬ 
electric  constants  correctly.  In  fact  the  dielectric  constants  are  overestimated 
by  the  first  parameter  set  (16.2  and  4.0)  and  are  underestimated  by  the  sec¬ 
ond  set  (8.1  and  2.9)  as  compared  to  the  experimental  values  of  11.0  and 
3.5.  This  discrepancy  may  be  due  to  the  different  description  of  cationic 
polarizability.  Bush’s  set  considers  rigid  Y  and  polarizable  Al  and,  in  con¬ 
trast,  Schuh’s  set  contains  rigid  Al  and  polarizable  Y.  On  the  basis  of  ionic 
radii,  Y  ion  is  expected  to  be  more  polarizable  than  Al  one.  As  it  is  shown 
from  our  calculations,  an  inclusion  of  both  Y  and  Al  polarizabilities  into  the 
model  allows  us  to  describe  dielectric  properties  of  YAG  very  well.  We  note 
here  that  the  accurate  representation  of  the  dielectric  constants  by  the  shell 
model  is  essential  in  getting  reliable  defect  energies  in  ionic  materials  such 
as  oxides. 

To  further  check  the  adequacy  of  the  derived  potential  parameters,  we 
examine  the  stability  of  the  YAG  lattice  with  respect  to  its  component  oxides: 


3Y2O3  T  5A/2O3  — 

>  2YzAl$Oi2i 

(6) 

A/2O3  +  3YAlOz  - 

Y^AhO\2- 

(7) 

The  reaction  enthalpies  come  out  to  be  -2.4  eV  and  -1.6  eV,  respectively,  per 
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YAG  formula  unit  indicating  that  YAG  is  indeed  a  stable  compound.  These 
values  could  be  compared  with  the  corresponding  reaction  enthalpies  for  YIG 
that  are  -1.6  eV  and  -2.9  eV,  respectively.21 

It  is  to  be  noted  here  that  the  structure  properties  of  orthorhombic 
perovskite-like  YAIO3  and  the  garnet  Y4AI4O12  compound  are  well  repro¬ 
duced  by  our  potential  model.  For  example,  the  calculated  lattice  constants 
of  YAIO3  (a=5.104  A,  b=5.260  A,  c=7.648  A)  are  in  good  agreement  with 
the  experimental  data  (5.179,  5.329,  7.370,  respectively42-43).  For  the  cubic 
Y4AI4O12  the  calculated  lattice  constant  of  12.106  A  is  only  0.9  %  larger  than 
experimental  value  of  11.989  A.44 


B.  Intrinsic  Defect  Structure 


For  calculations  of  native  point  defects  in  YAG,  the  lattice  was  simulated 
by  a  large  cluster  which  contained  more  than  250  ions.  The  ions  surrounding 
the  defect  in  the  cluster  were  allowed  to  relax  until  a  minimum  of  the  total 
energy  was  achieved.  The  rest  of  the  lattice  was  considered  in  terms  of  Mott- 
Litlleton  approximation.24  The  calculations  were  performed  using  the  GULP 
code.39 

A  vacancy  formation  was  simulated  as  a  removal  of  a  regular  lattice  site 
ion  to  infinity.  In  order  to  model  interstitials  one  ion  was  added  to  the  crystal 
from  infinity  and  was  located  in  the  empty  interstitial  position  in  the  lattice. 
In  the  case  of  oxygen  interstitials  we  probed  all  the  possible  positions  in  YAG 
which  are  referred  to  as  /  and  g.  For  cations,  the  b  positions  were  examined. 
The  calculated  formation  energies  for  the  isolated  basic  defects  (which  are 
collected  in  Appendix-Table  A. 2)  can  now  be  used  to  analyze  the  energetics 
of  Schottky  and  Frenkel  disorder  in  the  lattice.  For  reader’s  convenience 


9 


earlier  potential  sets17,40  and  the  defect  energies  borrowed  from17  as  well 
as  obtained  here  using  the  corresponding  set  of  pair  potentials  are  given  in 
Appendix  for  comparison. 

Since  YAG  is  a  complex  oxide  and  it’s  lattice  consists  of  three  sublattices, 
the  possible  Schottky-like  structure  disorder  can  be  written  in  the  form: 

Y AG  — >  I3A/5O12  +  3 Vy  +  5V^  +  12V0,  (8) 

Y AG  ->  A/203  +  2VZ  +  3Vq,  (9) 

YAG  — ■»  Y2O3  +  2Vy  +3 Vq,  QO) 


YAG  -»  YA103  +  Vy  +  V^l  +  3Vq,  (XI) 

where  we  use  the  notations  of  Kroger  and  Vink45  for  defects  and  the  symbol 
”YAG”  denotes  the  bulk  crystal. 

Here,  the  reaction  (8)  corresponds  to  the  true  Schottky  disorder  and 
remaining  ones  (9)-(ll)  are  given  as  examples  of  a  change  in  the  chemical 
composition  of  the  material.  They  can  be  interpreted  as  a  possibility  of  self¬ 
segregation  in  stoichiometric  YAG  crystal  by  the  Schottky-like  disorder  in 
either  of  the  cation  sublattices. 

The  Frenkel  disorder  can  be  written  as 


YAG 


+  AY', 


(12) 


YAG 


V;  +  Yr , 


(13) 


YAG  ->  Vq  +  O'i  (14) 

The  computed  formation  energies  (which  are  collected  in  Appendix- Table 
A. 3)  suggest  that  the  Schottky-like  disorder  appears  to  be  more  favorable 
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than  the  Frenkel  disorder.  Nevertheless,  even  Schottky  defect  formation  de¬ 
mands  rather  high  energy  (of  about  3. 9-4. 7  eV  per  defect).  Therefore,  we 
can  suppose  that  both  Frenkel  and  Schottky  defects  are  not  likely  to  occur 
in  YAG.  This  conclusion  is  in  agreement  with  the  results  has  been  reported 
by  Donnerberg  and  Catlow21  for  Yttrium  Iron  Garnet  (YIG)  crystals.  Fur¬ 
thermore,  large  enthalpies  of  reactions  (9)-(ll)  preclude  the  occurrence  of 
self-segregation  in  stoichiometric  YAG  crystal  as  also  observed  experimen¬ 
tally.46 

Next  we  consider  antisite  disorder  in  the  cation  sublattice  of  YAG  which 
can  be  described  as  follows: 


Yf  +  V(«)'  YX,m  +  ATy, 

(15) 

Yy  +  Al(d)‘  y;w  +  Al*y. 

(16) 

The  corresponding  antisite-pair  formation  energy  is  estimated  to  be  about 
0.9  eV  and  1.8  eV  per  defect  for  the  reactions  (15)  and  (16),  respectively, 
suggesting  that  the  antisite  disorder  is  preferred  over  Schottky  and  Frenkel 
disorder  in  YAG.  We  note  here  that  Y  prefers  to  be  located  at  the  Alfa)  site 
over  the  A1  (d)  site.  This  is  expected  on  the  basis  of  the  coordination  numbers 
of  these  ions  in  the  lattice:  Y  has  coordination  number  of  8  whereas  A1  (a) 
and  A1  (d)  ions  have  6  and  4  near-neighbor  oxygen  ions,  respectively  (see  Fig. 
1).  It  is  worth  to  note  here  that  antisite  pair  formation  in  perovskite-like 
compound  YAIO3  also  demands  very  low  energy  which  is  0.62  eV  per  defect 
indicating  that  these  defects  should  play  an  important  role  in  the  material. 

Experimental  studies  performed  by  Ashurov  et  a/.19  have  in  fact  suggested 
that  Y  can  occupy  its  regular  dodecahedral  sites  as  well  as  the  Al(a)  site  in 
YAG.  They  investigated  non-stoichiometric  Y3(Y;rAl2-i)Al30i2  crystals  by 
measuring  the  spectra  of  rare-earth  impurity  ions  in  the  lattice  and  have 
concluded  that  the  substitution  of  Y  at  the  Al(a)  site  introduces  the  lattice 
distortion.  Recent  XANES  and  EXAFS  measurements  by  Landron  et  al.20 
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have  investigated  the  local  order  around  yttrium  regular  and  antisite  loca¬ 
tions  in  Y2O3  rich  YAG.  As  it  was  shown,  an  increase  of  the  electronegativity 
of  the  site  ions  after  the  substitution  of  A1  cations  by  Y  implies  that  the  oxy¬ 
gen  bond  in  the  lattice  becomes  less  ionic.  The  results  in  shortening  of  the 
distance  between  Y  and  0  from  d(Y-0)=2.43  A  to  d(Y-O)=2.10  A  and  is 
now  close  to  the  mean  distance  d(Y-O)  of  yttrium  in  octahedral  position  in 
yttria  (Y2O3).  From  our  calculation,  the  substitution  of  Y  at  the  A1  (a)  site 
in  YAG  lattice  is  accompanied  by  the  decrease  of  the  distance  d(Y-O)  to  2.08 
A  which  is  in  excellent  agreement  with  the  experimental  value  of  2.10  A.20 


C.  Deviations  from  Stoichiometry 


Although  experimental  studies  have  long  ago  identified  the  presence  of 
non-stoichiometric  phases  in  YAG,  the  mechanism  by  which  an  excess  of 
yttria  or  alumina  can  be  accommodated  in  the  lattice  is  not  completely  un¬ 
derstood.  In  the  early  work  on  gallium  and  aluminum  garnets,36  excess  of 
yttrium  oxide  was  found  to  be  soluble  in  the  garnet.  Geller35  discussed  a  pos¬ 
sible  mechanism  of  the  solid  solution  in  the  garnet  structure  which  may  in¬ 
volve  vacancy  or  vacancy-interstitial  complex  instead  of  substitutional  atoms. 
Neiman  et  a/.14  reported  that  the  YAG  structure  can  exist  with  some  deficit 
of  AI2O3.  Several  studies  (see,  for  example,47,48)  have  also  identified  AI2O3 
inclusions  in  YAG  ceramics.  In  this  work,  we  make  an  aim  to  find  the  most 
probable  mechanism  of  accommodating  the  deviations  from  stoichiometry  in 
YAG  structure. 

For  Y2O3  excess  in  YAG,  we  can  write  the  following  equations  involving 
the  formation  of  vacancies,  interstitials  and  antisites: 

Y203  +  YAG->  2  Y^  +  AhOs  E  =  0.5eV,  (17) 
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Y203  +  YAG  -► 

2  Yr  +  30" 

E  =  36.2eV, 

(18) 

Y203  +  YAG  -* 

2YA103  +  2VZ  +  WS 

E  =  23.3eV, 

(19) 

Y203  +  YAG  -► 

lYAG  +  YZ  +  lVZ  +  VS 

E  =  7.4eV, 

(20) 

Y203  +  YAG  -* 

iYAG+irz 

E  =  -0.3eV, 

(21) 

Y203  +  YAG  -> 

AG  +  fir  +  |vj 

E  =  9.4eV, 

(22) 

Y203  +  YAG 

lYAG+i§VZ-+5Vo 

E  =  38.1  eld 

(23) 

The  subscript  A1  refers  to  the  regular  Al(0)  lattice  site  in  YAG  and  E  is  the 
calculated  enthalpy  of  the  given  reaction. 

Quite  similarly  we  can  describe  the  accommodation  of  excess  of  AI2O3: 


Al203  +  YAG  -> 

2  Aly  +  Y2O3 

E  =  2.9eY, 

(24) 

A/2O3  +  YAG  -> 

2  Air  +  30" 

£  =  52.5eV, 

(25) 

A/2O3  +  YAG  -»• 

2YA/O3  +  2  Vy  +  3  Vo 

E  =  18.9eY, 

(26) 

A/2O3  +  YAG  -* 

\YAG  +  Aly  +  \Vy  +  Vq 

£  =  5.7eV, 

(27) 

A/2O3  +  YAG  -> 

\Y  AG  +  f  A/y 

E  =  0.2eV, 

(28) 

A/2O3  +  YAG 

4  Y  AG  +  f  A/f  +  f  Yy" 

£  =  7.6eY, 

(29) 

A/2O3  +  YAG  -> 

fYAG+|Y;"  +  |Y0- 

E  =  10.7  eV. 

(30) 

Inspecting  the  list  of  the  reactions  and  energies  we  notice  that  only  re¬ 
action  (21)  has  negative  energy  and  is  therefore  exothermic.  It  leads  to  the 

conclusion  that  the  antisite  substitution  of  is  the  most  favorable  of  all 

the  possible  mechanisms  of  accommodation  of  extra  Y  ions  in  YAG  struc¬ 
ture.  This  energy  can  be  compared  with  the  energy  of  the  reaction  (17) 
which  describes  antisite  Y formation  accompanied  with  segregation  of 
AI2O3  and  shows  a  very  small  (endothermic)  reaction  enthalpy  indicating  a 
high  probability  of  its  occurrence.  Alumina  inclusion  formations47,48  may 
therefore  be  stimulated  by  yttria  surplus  in  YAG.  A  deficit  of  alumina  in  the 
garnet  structure,  that  was  postulated  by  Neiman  et  a/.,14  indeed  can  exist 
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in  YAG.  Comparing  the  energies  of  reactions  (21)  and  (28),  we  find  that  the 
Y^/(a)  antisite  formation  in  Y2O3  rich  YAG  is  slightly  more  preferable  over 
the  formation  of  Alf-  in  AI2O3  rich  YAG  although  the  latter  is  still  possible 
due  to  a  small  enthalpy  of  the  process.  Besides,  these  reactions  indicate  that 
extra  Y  ions  can  be  accommodated  in  YAG  lattice  easier  than  A1  ions.  This 
fact  is  known  from  experiments.  For  example,  an  excess  of  Y2O3  up  to  2% 
was  reported  while  an  surplus  of  AI2O3  in  YAG  was  seen  only  up  to  0.5% 
(see,  for  instance,18). 

The  formation  of  defect  complexes  like  an  interstitial  and  the  correspond¬ 
ing  vacancy  (see  reactions  (22)  and  (29))  requires  much  more  energy  with 
respect  to  the  antisite  formation  (reactions  (21)  and  (28)).  The  required  en¬ 
ergy  for  the  formation  of  interstitial  defects  described  by  reactions  (18)  and 
(25)  are  also  very  high  as  it  was  expected  from  Frenkel  disorder  energetics 
(see  Section  III.B).  The  segregation  in  the  perovskite-like  phase  (reactions 
(19)  and  (26))  in  YAG  is  also  expected  to  have  a  low  probability  due  to  very 
high  energy  of  the  reaction. 

The  reaction  (20)  represents  a  mechanism  of  the  defect  accommodation 
predicted  by  Neiman  et  al.  in14  while  equation  (27)  is  composed  by  analogy. 
And  the  reactions  (23)  and  (30)  describe  the  process  suggested  by  Sakaguchi 
et  al.  in.18  From  our  calculations,  both  of  these  mechanisms  are  not  pre¬ 
dicted  to  take  place  in  YAG  crystal.  The  calculated  results  imply  that  the 
defect  equilibrium  in  an  excess  material  has  different  nature  than  that  con¬ 
sidered  in  these  studies.  Excess  of  neither  of  Y2O3  nor  AI2O3  appears  to 
be  source  of  oxygen  vacancies  in  pure  YAG.  It  should  be  noted  here,  that 
an  interaction  between  defects,  that  may  be  important  under  the  conditions 
of  large  deviations  from  stoichiometry  in  materials,  was  not  included  in  the 
present  study. 
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IV.  CONCLUSIONS 


In  the  present  theoretical  study,  atomistic  calculations  of  a  formation  of 
native  point  defects  in  Yttrium  Aluminum  Garnet  crystal  are  performed. 
The  Schottky  and  Frenkel  defect  formation  energies  are  calculated  as  well  as 
antisite  substitutions  in  the  cation  sublattice.  The  calculated  results  show 
that: 

1.  intrinsic  disorder  in  stoichiometric  YAG  is  dominated  by  antisites; 

2.  antisite  substitution  Y xAl^  causes  a  distortion  in  the  crystalline  lattice 
shortening  the  Y-0  bond  length  whose  calculated  value  is  in  excellent  agree¬ 
ment  with  the  EXAFS  measurements; 

3.  excess  of  yttrium  oxide  Y2O3  is  energetically  more  favorable  with  respect 
to  aluminum  oxide  AI2O3  surplus  in  YAG; 

4.  formation  of  antisites  in  non-stoichiometric  YAG  is  the  most  probable 
mechanism  of  accommodating  the  deviations  from  stoichiometry. 
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TABLE  1.  Interatomic  potential  parameters  and  shell  constants.  A,  p  and  C  are  parameters 
in  the  equation  for  the  Buckingham  potential  (see  in  the  text,  eq.  (4));  Y+  is  charge  of  the 
core;  K+  is  core-shell  spring  constant. 


Ion 

A  (eV) 

P  (A) 

Y+(|e|) 

K+  (eVA'2) 

C  (eVA'6) 

o2  -  o2- a) 

22764.000 

0.149 

0.8481 

74.92 

27.88 

y3+-  o2 

2036.8379 

0.3103 

-0.251 

46.7288 

0.0 

ai3+  -  O2- 

741.9007 

0.3566 

0.043 

40.8618 

0.0 

a)  potentials  for  O2'-  O2-  interaction  are  taken  from  G.V.Lewis  and  C.R.A.Catlow  (Ref.  41). 


TABLE  II.  Structural  data  for  YAG  calculated  using  interatomic  potential  parameter  sets  of: 
l.the  present  study;  2.  Schuh  et  al.  17  ;  3.  Bush  et  al.  40. 
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APPENDIX 


Table  A.l  shows  the  interionic  potential  and  shell  model  parameter  sets 
of  Schuh  et  al.17  and  Bush  et  al.40  Table  A. 2  gives  basic  isolated  defect 
formation  energies  obtained  using  these  parameter  sets  along  with  our  derived 
set. 

One  of  the  most  important  advantages  of  the  shell  model  is  a  correct 
description  of  optical  modes  of  vibration,  and  as  a  result,  a  value  of  the  high 
frequency  dielectric  constant.  As  it  was  shown  (see  sec.  III.B),  this  value 
may  be  sensitive  to  the  treatment  of  polarization  of  cations  in  the  model. 

Now  we  can  analyze  an  influence  of  polarization  on  defect  formation  en¬ 
ergy.  Energies  of  Schottky-like  and  Frenkel  disorder  in  YAG  structure  are 
collected  in  Table  A. 3  from  which  it  is  seen  that  the  choice  of  potential  does 
not  change  the  main  conclusion  about  preference  of  the  disorder  in  YAG. 
However,  a  more  careful  look  on  the  energetics  of  reactions  reveals  some  sort 
of  contradiction.  The  formation  of  Frenkel  Y  defects  (reaction  (13))  demands 
the  highest  energy  using  Schuh’s  set,  whereas  using  parameters  of  Bush  the 
highest  energy  comes  out  to  be  for  Frenkel  Al(d)  defects  (reaction  (12)).  The 
lowest  energy  among  Frenkel  pairs  corresponds  to  Al  (d)  pair  formation  in 
the  Schuh’s  set  quite  contrary  to  the  Bush’s  results. 

Furthermore,  the  probability  of  segregation  of  alumina  and  yttria  (reac¬ 
tions  (9)  and  (10))  is  found  to  be  comparable  by  Schuh’s  set.  But  Bush’s 
set  predicts  a  large  difference  between  these  reactions.  Also,  from  Bush’s 
parameters,  segregation  of  yttrium  oxide  (reaction  (10))  is  only  0.4  eV  less 
favorable  than  antisite  pair  formation  in  contrast  to  the  Schuh’s  results. 

These  results  lead  to  the  conclusion  that  defect  formation  energy  is  very 
sensitive  to  the  treatment  of  polarization  effects  in  the  lattice.  Therefore,  in 
order  to  get  accurate  results  on  energetics  of  the  defect  formation  in  YAG 
one  should  take  carefully  into  account  polarization  of  ions  in  the  lattice. 
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TABLE  A.I  Interatomic  potential  parameters  and  shell  constants  of  Schuh  et  al.  17  ; 

A,  p  and  C  are  parameters  in  the  equation  for  the  Buckingham  potential;  Y+  is  charge  of  the 
core;  K+  is  core-shell  spring  constant;  C  is  attractive  contribution. 


Ion 

A  (eV) 

P  (A) 

Y+(|e|) 

K+  (eVA'2) 

C  (eVA6) 

O2'-  o2- 

22764.000 

0.149 

0.8481 

74.92 

27.88 

Y3+-  O2' 

1345.1 

0.3491 

-0.251 

71.70 

0.0 

ai3+-  o2- 

1469.3 

0.2991 

3.00 

HBH 

0.0 

Interatomic  potential  parameters  and  shell  constants  of  Bush  et  al.  40  . 


Ion 

A  (eV) 

P  (A) 

Y+(|e|) 

K+  (eVA‘2) 

C  (eVA'6) 

r!» 

O 

1 

r-» 

O 

25.41 

0.6937 

0.513 

20.53 

32.32 

Y3+-  O2- 

1519.279 

0.3291 

3.00 

hbh 

0.0 

> 

U) 

+ 

l 

o 

2409.505 

0.2649 

0.043 

403.98 

0.0 

TABLE  A.2.  Intrinsic  defects  in  YAG.  Defect  fiormation  energies  obtained  using  different 
parameter  sets. 


Vacancy  formation  energy  (eV) 

Parameters  set 
of 

0 

Y 

Al  (a-site) 

Al  (d-site) 

present  study 

21.54 

49.42 

53.76 

53.44 

Schuh  et  al. 

22.35a> 

46.5a) 

60.24a) 

58.9a) 

Bush  et  al. 

17.63 

47.74 

59.94 

65.10 

Interstitial  formation  energy  (eV) 

0 

O 

Y 

Al 

Parameters  set 

Symmetry  position 

of 

8- 

(1/8,  x,  1/4-x) 

f: 

(x,  0,  1/4) 

b: 

(1/8,  1/8,  1/8) 

b: 

(1/8,  1/8,  1/8) 

present  study 

-11.38 

-11.79 

-36.85 

-40.72 

Schuh  et  al. 

-11.35 

-11.06 

-25.62a) 

-45 . 1 a) 

Bush  et  al. 

-8.33 

-7.84 

-33.51 

-47.48 

Antisite  formation  energy  (eV) 

Parameters  set  of 

Al  in  Y  site 

Y  in  Al(a)  site 

Y  in  Al(d)  site 

present  study 

-0.69 

2.38 

4.16 

Schuh  et  al. 

-10.45 

13.53 

14.29 

Bush  et  al. 

-8.93 

11.52 

12.78 

3l\ 


a)  from  L.Schuh  et  al.  17 


TABLE  A.3.  Reaction  energies  (per  defect)  in  eV  for  Schottky-like,  Frenkel  and  antisite 
disorder  reactions  obtained  using  different  parameter  sets.  For  details,  see  the  text. 


Reaction 


Parameters  set  of 

Present  study 

Schuh  et  al. 

Bush  et  al. 

Schottky-like 


(8) 

4.26 

(9) 

4.70 

(10) 

3.88 

(ID 

4.21 

5.32a) 


4.93 


4.91 


Frenkel 

(12):  Al  from  (a)  site 

6.65 

7.57  a) 

6.23 

(12):  Al  from  (d)  site 

6.36 

6.88  a) 

8.81 

(13):  Y 

6.28 

10.44  a) 

7.12 

(14):  0  with  (g)  interstitial 

5.08 

5.50 

4.65 

(14):  O  with  (f)  interstitial 

4.88 

5.65 

4.90 

Antisite  pair 

(15) 

0.9 

1.54 

1.30 

(16) 

1.8 

1.92 

1.95 

the  value  is  taken  from  Schuh  et  al. 


Figure  caption  for 

Atomistic  modeling  of  defects  in  Yttrium  Aluminum  Garnet 

crystals 

by  M.M.Kuklja  and  R.Pandey. 

FIG.  1.  Schematic  representation  of  the  arrangement  of  YAG  in  terms  of 
the  garnet  structure.  A  fragment  of  a  cluster  is  centered  at:  a),  yttrium  atom 
which  is  dodecahedrally  coordinated  (Y(c));  b).  aluminum  atom  which  is 
octahedrally  coordinated  (Al(a));  c).  aluminum  atom  which  is  tetrahedrally 
coordinated  (Al(d))  in  the  crystal. 
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Fig, 
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